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We present a formalism for studying the influence of dispersive and absorbing 
dielectric bodies on a radiating atom in the framework of quantization of the phe- 
nomenological Maxwell equations for given complex permittivities of the bodies. 
In Markov approximation, the rate of spontaneous decay and the line shift asso- 
ciated with it can then be related to the complex permittivities and geometries of 
the bodies via the dyadic Green function of the classical boundary value problem 
of electrodynamics - a result which is in agreement with second-order calcula- 
tions for microscopic model systems. The theory is applied to an atom near a 
planar interface as well as to an atom in a spherical cavity. The latter, also known 
as the real-cavity model for spontaneous decay of an excited atom embedded in 
a dielectric, is compared with the virtual-cavity model. Connections with other 
approaches are mentioned and the results are compared. 



1 Introduction 

It is well known that the quantum statistics of photons and atoms can essentially be 
influenced by the presence of macroscopic bodies. In principle, such bodies can be 
included in the quantum mechanical description of the system as a part of matter to 
which the radiation field is coupled and treated microscopically, using, e.g., a canonical 
Hamiltonian formalism and perturbation theory (mostly up to second order). Apart 
from the fact that in practice the involved calculations can only be performed for rel- 
atively simple systems, each system requires its own calculation in general, because of 
the specific microscopic properties that must be taken into account. 

In classical optics the effect of macroscopic bodies is commonly described phe- 
nomenologically in terms of a spatially varying permittivity e(r,ui) [or refractive in- 
dex n(r, u/) = yj e(r, ui)] which is a complex function of frequency in general. Recently a 
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scheme for transferring this concept to quantum theory has been developed [1,2,3] which 
is consistent with standard QED as well as with the dissipation-fluctuation theorem. 
The advantage of such an approach is that it is valid for arbitrary linear dielectrics. 
Knowledge is needed only about the dispersion-absorption profile of the medium ex- 
pressed in terms of the (measurable) permittivity, from which the dyadic Green function 
of the classical problem can be calculated. The electric and magnetic fields can then 
be expressed in terms of the classical Green function and a continuum of fundamental 
bosonic fields. 

In the present paper we use that quantization scheme for studying the problem 
of a radiating atom in the presence of dispersive and absorbing dielectric matter. In 
particular, we consider the effect on the spontaneous decay of an excited atom and 
calculate the decay rate and the line shift. With regard to potential applications in 
scanning near-field optical microscopy, we present results for an atom near a planar 
interface. Further, we derive a correct generalization for absorbing media of the real- 
cavity model and the virtual-cavity model of spontaneous decay in dielectrics. 

The paper is organized as follows. In Section 2 the quantization scheme for the 
electromagnetic field in linear Kramers-Kronig dielectrics is briefly reviewed, and general 
formulas for the decay rate and the line shift of an atom in the vicinity of a dielectric 
body are derived. In Section 3 a planar interface is considered, and Section 4 is devoted 
to a spherical cavity, which may serve as a model for spontaneous decay in media. The 
corresponding virtual-cavity model is considered in Section 5, and a comparison with 
the real-cavity model is given. Finally, a summary is given in Section 6. 



2 Basic equations 

Our analysis is based on the phenomenological quantization of the electromagnetic field 
in linear Kramers-Kronig dielectrics as developed in [1,2,3]. Without external sources, 
Maxwell's equations in (temporal) Fourier space read 

V-B(r,w) = ) (1) 

V- [e e(r,w)E(r,w)] =£(r,w), (2) 

V x E(r, lo) = iuV x B(r, w), (3) 

V x B(r, w) = -i^e(r, w)E(r, w) + /i j(r, w), (4) 

cr 

with e(r, u>) = en(r,u>) +iei(r,u>) being the spatially varying, complex permittivity of 
the medium which satisfies the Kramers-Kronig relations. In Eqs. (2) and (4), p(r,u) 

and j (r, u) , respectively, are the operator noise charge and current density operators 
that are associated with absorption according to the dissipation-fluctuation theorem. 
They are related to the noise polarization P JV (r, u) as 

j(r,w) = -iuP N {r,u), (5) 
p(r,u>) = -V-F(r, W ), (6) 
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and hence satisfy the equation of continuity. By introducing an infinite set of bosonic 
field operators f(r,u), all electromagnetic field operators can be expressed in terms of 
them using the fundamental relation 



j(r, u) = w\j — e/(r, w) f (r, w). 



(7) 



In particular, Maxwell's equations imply that the electric field operator may be written 
in terms of the classical dyadic Green function as 



E k (r,uj) = ifjQ J d 3 sG fcfc /(r,r',w)j fe ,(r',w), 
where G kk > (r, s, w) satisfies the partial differential equation 



dldl-8 ik [A r + —e{r,uj) 



G k i{r,v',w) = 5 u 8{v-v'). 



Note that the relation 



/ 



d 3 s -J ej(s, w)G k i (r, s, u)G%, l (r', s, w) = Im G kk > (r, r', w) 



(8) 



(9) 



(10) 



is valid, which can be checked by adding the complex conjugate of Eq. (9) to itself and 
integrating by parts. 

Integrating Eq. (8) with respect to to leads to the operator of the electric field and, via 
Maxwell's equations, to the magnetic field operator, with the equal-time commutation 
relations 



£fc(r),B;(r'; 



-ihe k[m d r m 5(r - r'), 



E k (T),Ei(r f ) = S fe (r),A(r') 



= 0, 



(11) 
(12) 



which has been shown [3] to be valid for arbitrary linear dielectrics that are consistent 
with the Kramers-Kronig relations. 

If there is an additional (two-level) atom present, then the total Hamiltonian of the 
system in dipole and rotating wave approximation can be given by 

d 3 r dw?iwf t (r,a;)-f(r, u>) + V] Tiuj a A aa - iuj 2 iA 21 A (+) (r A ) d 2 i + H.c. 
J o a= i L 

(13) 

where we have chosen the Weyl gauge for the vector potential, and hence the scalar 
potential is set equal to zero. It can then be shown that in Markov approximation the 
atomic (flip) operators A aa > = \a)(a'\ satisfy the following Langevin-type differential 
equations [4]: 



122 



-TA 



22 



4yAW(r,,i).d 21 |H,. 



(14) 

An = -A 22 (15) 
i 2 i = [i(w 21 -M-|r]i 21 + ^A^(r A ,i) -dai^ -A n ), (16) 
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where the rate of spontaneous decay of the excited atomic state |2) is given by 



r = ^ImG t ,(r Al r A ,, A ), (17) 



and the Lamb shift contribution reads 

Im G k i(r A ,r A , u) 



he c 2 



ReG k i(rA,r A , wa) - - / doj Im Gkl 
7T Jo w 



+ WA 



(18) 



[/ife = (c?2i)fe, wa =W2i]- Note that the coincidence limit of the real part of the vacuum 
Green function is infinite and regularization is required. Here and in the following the 
vacuum Lamb shift is thought of as being included in the atomic transition frequency so 
that Sui only results from the medium. With respect to Green functions of the type given 
in Eq. (21), with G\i(ya, ?a, wa) being the vacuum Green function, Gki(^A, ?A, wa) in 
Eq. (18) can therefore be replaced with R k i{rA, ^A, ^a), and Eq. (17) reduces to 

r = r + 2 "^f l lmR k i{v A ,VA,0JA), (19) 

where r = uj a ^ 2 / (3nheoc 3 ) is the spontaneous emission rate of the atom in free space. 

Equations (17) and (18) hold for all linear dielectrics and arbitrary geometries. Note, 
that Eq. (17) can be rewritten, on using the relation (10), as 



ImG t ,(r,r»^- W ') = l^H 01 l^ fe(r ' (r '' "'A |0) ' (20) 

which is nothing but a consequence of the consistent introduction of the noise polariza- 
tion according to the dissipation-fluctuation theorem. It should be pointed out that for 
an atom in vacuum the imaginary part of the Green function for arbitrary geometries 
of the surrounding dielectric is purely transverse and, most importantly, finite. 
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3 Planar interface 

Let us first consider the spontaneous decay of an excited atom near an absorbing planar 
dielectric surface. This configuration with real refractive index has been studied exten- 
sively in the literature in connection with Casimir and van der Waals forces (see [5] 
and references therein) as well as applications in scanning near-field optical microscopy 
(SNOM) (see [6] and references therein). To be more specific, we consider two infinite 
half-spaces (volumes Vi and V2) with a common interface such that e(r, uj) = 1 if r 6 Vi, 
and e(r, uS) = e(uj) if r G V2 and assume that the atom is placed in Vi at a distance z 
from the interface. The Green function in V\ for this configuration can be given by [2] 

G kk ,(v,v',uj) = G k / k ,(r,r',ij) + R kk ,(v,v',uj) (r.r'eVi), (21) 

where G\[ k , (r, r', u>) is the vacuum Green function, and R kk > (r, r', ui) is the reflection 
Green function. As already mentioned, for calculating the decay rate and the line shift 
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we only need the reflection Green function with both spatial arguments at the position 
of the atom, R/-k' = Rkk' (ya,ya,u). Using the formulas given in [2], we obtain 

R,, = R„„ = -5— o / dkkpe^r^k) + — dk — — r s (fc), (22) 
»7T3 Jo on J Q [3 



2i(3z 



(q = u)/c, (3= \J q 2 — k 2 ). The functions r p (k) and r s (k) are the usual Fresnel reflection 
coefficients for p- and s-polarized waves. The integrals can be evaluated asymptotically 
in the case when qz <C 1, i.e. when the distance of the atom from the surface is smaller 
than the transition wavelength, which is in agreement with the Markov approximation. 
Expanding in powers of z leads to [e = e(w), n = y 7 e(u>)] 

1 n 2 - 1 1 (n-1) 2 iq (n-l)(2n-l) 
zz 16nq 2 z 3 n 2 + 1 8nz n(n + 1) 12tt n(n+l) k j ' V ; 



2 l07rz + 1 37rn + 1 



Inserting Eqs. (24) and (25) into Eq. (19) yields 

The term proportional to z~ 3 exactly agrees with the formula for the leading term 
derived in [7] in the framework of a microscopic approach, after lengthy and involved 
diagonalization of the model Hamiltonian used and computation of perturbation sc- 
ries up to second order. It is worth noting that the phenomcnologically introduced 
permittivity and the Green function that is based on it already contain the relevant in- 
formation, so that the effort is reduced drastically. Needless to say that the microscopic 
treatment given in [7] is far from a genuine ab initio calculation. 

SNOM detects surface corrugation or impurities via the changes in the line shift 
and the line width of an atomic transition far from medium resonances [6]. The high 
sensitivity of the method results from the cubic dependence of the line shift on the 
inverse distance of the probe atom from the surface. Disregarding absorption the line 
width is weakly distance-dependent [cf. Eqs. (19) and (24)], and the strong dependence 
on distance of the line shift is essentially determined by the first term in Eq. (18), 

The situation changes drastically in an absorption band or in the vicinity of it. From 
Eq. (26) it is seen, that in an absorbing regime the line width also becomes proportional 
to the inverse cube of the distance. Moreover, the line shift can substantially increase. 

To demonstrate the possibility of improving the resolution of SNOM near a medium 
resonance, we have calculated numerically the spontaneous decay rate as well as the 
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line shift using the single-resonance model permittivity 

(0.46^ r ) 2 
e{LO) = 1 + — j 5 : — . 



(28) 



From Figures 1 and 2 one observes that near the medium resonance both the decay 



120 
100 
80 
60 
40- 
20 

ot. 



0.2 0.4 



0.6 0.8 



1 1.2 




Fig. 1. Decay rate T/To (/i = fi z ) for = 0.05 = 1.29] as a function of frequency 

(left figure) for 2tvz/\a = 0.1 (full line) and 2ttz/\a =0.3 (dashed line), and as a function of 
distance (right figure) for lua/^t = 1 (full line) and lja/^t = 0.5 (dashed line). 
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Fig. 2. Line shift Slo/Tq (h = Hz) for ^/uja 
(left figure) for 2itz/\a =0.1 (full line) and 2ttz/\a =0.3 (dashed line), and as a function of 
distance (right figure) for lua/ut = 1 (full line) and loa/^t = 0.5 (dashed line). 
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rate and the line shift exceed substantially the values realized in standard SNOM. 
Moreover, the resolution of scanning is improved, as can be seen from the vertical 



line- width resolution 



ST 



1 



oc^ , ;^ A \ {2 (29) 
Sz z 4 \e(oj A ) + 1| 2 

(a similar formula holds for the vertical line-shift resolution). Hence we are left with 
the proposal of operating in an absorbing regime near a medium resonance in order to 
improve the sensitivity of measurement. 



4 An excited atom in a spherical cavity 

Another interesting geometry is the spherical one. In particular, an atom in an empty 
sphere Vr of radius R inside an otherwise continuous dielectric medium can serve, in 
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a sense, as a model for spontaneous decay in media (real-cavity model). For r,r' G V R 
the Green function can again be written in the form of Eq. (21), where the reflection 
tensor reads in dyadic notation as [8] 



i?(r,r» 



2n + 1 (n — to)! 
4ttc ^— ' ^— : ' ^— ' 1 n(n + 1) (n + to)! 

e,o n— 1 m— ^ v 7 v 7 



EE 



E 



(2 — cWi) 



Cf (u;)M >m (r, -^)M >m (r', ^) +C^(u)N e nm (r, >m (r' , ^ 



(30) 



Here, A / f = „ m (r, fc) and AT = „ m (r, fc) are the (even and odd) vector Debye potentials, and 
the quantities C^f ,N (u>) are the generalized reflection coefficients. Combining Eqs. (19) 
and (30) (for r = r' = r A ), after some algebra we derive 



where 



r= [1 + ReCf (wa)] r °> 



[i + z(n + 1) - zz 2 n - z 3 n 2 /(n + 1)] e* 



z(cosz + in sin z) + iz 2 ncosz — z 3 (cosz — msinz)n 2 /(n 2 — 1) 



(31) 



(32) 



(z = Ruj/c; for details, see [4]). Examples are shown in Figs. 3 and 4. Note that the 
decay rate is solely related to the transverse electromagnetic field, r = T 1 - . 
Restricting our attention to small cavities, loaR/c-^ 1, expansion yields 



r = 



\\2e(ui A ) + 1\ 2 \lo a R 
9rj(uj A ) 



+ 



9ej(u; A ) [28| £ (^)| 2 + 12e R (u; A ) + l] / c 



|2e(^) + l| 4 

9k(u>a)£i{ua) 
\2e{u A ) 



5\2e(Lu A ) + 1\ 4 
4|e(w j4 )| 4 + 4e R (uj A )\e(uj A )\ 2 + e R (co A ) - e 2 (uj A ) 



uj a R 



II 



4\e(LU A )\ 2 + 2e R (LU A ) 



>v a + o{R) 



(33) 



(e R = Kce, ej = Ime). When absorption is fully disregarded, i.e., Ime = 0, Eq. (33) 
reduces exactly to the well-known Glauber-Lewenstcin formula [9] 



3n 2 



2n 2 



1 



nr . 



(34) 



Needless to say that when setting e = 1 , then the free-space spontaneous emission rate is 
recovered. In absorbing media the decay rate sensitively depends on the cavity radius. 
Note that the terms ~ i?~ 3 and ~ R^ 1 result from absorption. In particular, the term 
<~ R~ 3 can be regarded as describing nonradiative energy transfer from the guest atom 
to the constituents of the dielectric via dipole-dipole interaction. It should be pointed 
out that the above derived rate formula for absorbing media strongly contradicts the 
rate formula suggested in [10]. 



8 



S. Scheel, L. Knoll, D.-G. Welsch 



5 Virtual-cavity model of spontaneous decay 

In the (Clausius-Mosotti) virtual-cavity model it is assumed that the field outside the 
small sphere is not modified by the sphere, and the local field inside is given by the 
(smeared) macroscopic field corrected by adding a contribution from the polarization 
of the sphere. In quantum theory, one has to be careful with the definition of the local 
field, because of the nonvanishing noise polarization even in the zero-temperature limit. 
In particular, simple multiplication of the macroscopic electric field with a correction 
factor as suggested in [10,11] leads to a wrong expression for the rate of spontaneous 
decay and misinterpretations concerning the action of transverse and longitudinal fields. 
Introducing the local field as E'(r,w) = E(r, u) +P(r, w)/(3eo), from Eqs. (2) and (4) 
together with Eqs. (5) - (7) it follows that [4,12] 



E'(r,o/) = \[e(v,u) + 2]E(r,w) + -Lp w (r, W ), 



with 



-j(r,w) 



hen 



e/(r,w)f(r,w). 



(35) 



(36) 



The spontaneous decay rate can now be calculated using Eqs. (17) and (20) to- 
gether with the local field defined in Eq. (35) and the Green function for bulk material, 
Gjjf/ (r, r', loa)- However, for r,r' -^>va the imaginary part of G^f k , (r, r', u>a) becomes 
singular, which reflects the fact that the dielectric is thought of as being smeared out 
even over the region where the atom is located. In order to overcome this difficulty, 
rcgularization is required. A small distance between the spatial arguments r and r' can 
be kept, and the decay rate can be obtained by averaging over an appropriately chosen 
sphere of radius R (virtual cavity), 



3heoC 2 



ImG£f,(r,r',^) + 



9he 



e I (aj A )6 kk >6(r - r') 



6/(0^) Re 



(37) 



(for details, see [4]). Obviously, different averaging procedures leads to (slightly) differ- 
ent numerical prefactors in the i?-dependent terms [4,12]. 

Let us consider a small cavity (\RuiAn/c\ <C 1). Averaging with respect to both 
spatial arguments separately over the sphere with equal weight, we obtain [12] 



r = r 1 - + r" 



(38) 



54 \loaR 



+ ei{uj A ) [ejjK) + 2] 



[A 




2k(uja) 






9 



£(wa) 



2ej(cu A ) 



>T Q + 0(R), 



(39) 
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Here, T 1 - and r", respectively, are related to the transverse and longitudinal parts of 
the electromagnetic fields. It should be pointed out that the appearance of r'l is due 
to the fact that (in contrast of the real-cavity model) the atom is not strictly located 
in vacuum and the imaginary part of the longitudinal Green function (at equal space 
points) does not vanish for absorbing media. 




0.8 0.9 1 1.1 1.2 0.8 0.9 " 1 " 1.1 1.2 



Fig. 3. Spontaneous decay rate F /Fq for R = 0.02Aa and 7 = 0.05u>t (left figure) and 7 = 
0.2u>t (right figure). The solid lines correspond to the real-cavity model, Eq. (33), and the 
dotted lines correspond to the virtual-cavity model, Eqs. (39) and (40), with the broken line 
indicating the transverse part (39) . 
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Fig. 4. Spontaneous decay rate F/Fo for R = 0.2Xa and 7 = 0.05cjt (left figure) and 7 = 0.2ljt 
(right figure). The solid lines correspond to the real-cavity model, Eq. (33), and the dotted 
lines correspond to the virtual-cavity model, Eqs. (39) and (40), with the broken line indicating 
the transverse part (39). 

Examples of T and T 1 - are shown in Figs. 3 and 4. The results differ from those 
in [10,11] in the contribution of the noise polarization P^, which is disregarded in 
[10,11]. It is worth noting that the effect of the noise polarization has recently been 
fully confirmed on the basis of a microscopic approach [13]. The effect vanishes for 
nonabsorbing media (ej = 0), and Eq. (39) yields the well-known rate formula 

^M 2 ^) 2 "^ (41) 
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6 Summary 

We have shown how quantization of the phcnomenological electrodynamics in linear 
media can be most easily employed to study radiating atoms in the presence of disper- 
sive and absorbing dielectric bodies and to calculate the line widths and line shifts of 
the atomic transitions. It is worth noting that the approach reproduces typical results 
that have been derived using quite different concepts including involved microscopic 
ones. The approach also shows that the rate formulas given in [10,11] for the spon- 
taneous decay of excited atoms in absorbing media are not correct, because they are 
based on wrong respective incomplete medium-assisted vacuum noise. In particular, 
we have considered a planar interface and a spherical cavity, which may serve as a 
model of spontaneous decay in media, and we have compared it with the corresponding 
(Clausius-Mosotti) virtual-cavity model. The results reveal that when an excited atom 
approaches an absorbing dielectric body, then the rate of spontaneous decay becomes 
proportional to the inverse cubic distance of the atom from dielectric body [see, e.g., 
Eqs. (26) and (33)]. In microscopic models, terms of that type have been shown to cor- 
respond to nonradiative energy transfer from the atom to the medium via dipole-dipolc 
interaction, i.e., a second-order process with no photon in the final state. With regard 
to applications, the results show that absorption offers the possibility of improving the 
resolution of SNOM. 
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